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Abstract 



We construct the metric of new multi-black hole configurations on a d-dimensional cylinder 
B^ -1 x S 1 , in the limit of small total mass (or equivalently in the limit of a large cylinder). 
These solutions are valid to first order in the total mass and describe configurations with 
several small black holes located at different points along the circle direction of the cylinder. 
We explain that a static configuration of black holes is required to be in equilibrium such that 
the external force on each black hole is zero, and we examine the resulting conditions. The 
first-order corrected thermodynamics of the solutions is obtained and a Newtonian interpreta- 
tion of it is given. We then study the consequences of the multi-black hole configurations for 
the phase structure of static Kaluza-Klein black holes and show that our new solutions imply 
continuous non-uniqueness in the phase diagram. The new multi-black hole configurations 
raise the question of existence of new non- uniform black strings. Finally, a further analysis of 
the three-black hole configuration suggests the possibility of a new class of static lumpy black 
holes in Kaluza-Klein space. 
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1 Introduction 

Black holes in four-dimensional General Relativity have a very simple phase structure. The 
uniqueness theorems for pure gravity assert that the only possible stationary black hole solu- 
tion for a given mass and angular momentum is the Kerr black hole. 

For higher-dimensional General Relativity, the situation is vastly different. In particular, 
if we imagine that we live in a world which is five dimensional with the extra dimension curled 
up on a circle, the relevant black hole solutions are those which asymptote to four-dimensional 
Minkowski-space times a circle (A4 x S 1 ), i.e. the five-dimensional Kaluza-Klein space-times. 
The phase structure of such black holes has been shown to be very rich and contains phases 
with event horizons of different topology and even phases where Kaluza-Klein bubbles are 
attached to black holes [H EJ [3] . More generally, we get a similarly rich phase structure for 
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the case of black holes asymptoting to (i-dimensional Minkowski-space times a circle (A4 d xS l ) 
with d > 40 The spatial part of this space-time is a d-dimensional cylinder x S 1 . 

The two static black hole phases which most obviously should appear for Ai d x S 1 are the 
localized black hole phase, which for small mass behaves as a d+ 1 dimensional Schwarzschild 
black hole, and the uniform black string corresponding to a d-dimensional Schwarzschild 
black hole times a circle. For the uniform string phase, the metric is known exactly. The 
most interesting feature of the uniform string is the Gregory-Laflamme instability [HJ [9] which 
is a long wave-length gravitational instability of the solution (see [3] for a review). From 
this instability, it follows that the uniform string for a certain mass has a marginal mode. 
From this marginal mode emanates a new branch of solutions which are non-uniform strings, 
i.e. solutions with same topology of the event horizons as the uniform strings but without 
translational symmetry around the circle. These new solutions have been studied numerically 

in pmini nana Hang. 

For the localized black hole phase, here dubbed the black hole on the cylinder phase, the 
metric is not known analytically. However, for small black holes on the cylinder the first 
order part of the metric has been found [TBI El OS IH] and also the second order solution 
has been studied [20], [21] . Finite-size black holes on the cylinder have instead been studied 
numerically [22\ [251 [24"] . This study has revealed the interesting result that the black hole 
on the cylinder phase meets the non-uniform string phase in a topology changing transition 
point [251 [261 [271 15]. 

In this paper, we find and study new solutions for multi-black hole configurations on the 
cylinder. These solutions describe configurations with several small black holes located at 
different points along the circle direction of the cylinder x S 1 . The location of each 

black hole are such that the total force on each of them is zero, ensuring that they are in 
equilibrium. It is moreover necessary for being in equilibrium that the black holes are all 
located in the same point in the M. d ~ l part of the cylinder. 

The metrics that we find are solutions to the Einstein equations to first order in the mass. 
More precisely, we work in a regime where the gravitational interaction between any one of 
the black holes and the others (and their images on the circle) is small. Thus, our solutions 
describe the small mass limit of these multi-black hole configurations on the cylinder, or 
equivalently they can be said to describe the situation where the black holes are far apart. 
The technique used for solving the Einstein equations is the one developed in [T7J for small 
black holes on the cylinder based on an ansatz for the metric found in |16| . 

A subset of the multi-black hole configurations have already been studied in the literature. 
These are the so-called copies of the black hole on the cylinder solutions |28[ [29l [T7] . This 
class of solutions corresponds to the special situation in which a number of black holes of the 
same size are spread with equal distance from each other on the circle. 

The existence of these new solutions have striking consequences for the phase structure 
of black hole solutions on Ai d x S 1 . It means that one can for example start from a solution 

1 Note that the case M 3 x S 1 studied in [J] [51 [7] is different due to the high amount of symmetry. 
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with two equal size black holes, placed oppositely to each other on the cylinder, and then 
continuously deform the solution to be arbitrarily close to a solution with only one black 
hole (the other black hole being arbitrarily small in comparison). Thus, we get a continuous 
span of classical static solutions for a given total mass. This means for static black hole 
solutions on A4 d x S 1 we have in fact a continuous non-uniqueness of solutions. Continuous 
non-uniqueness for black holes has also been found when one attaches Kaluza-Klein bubbles 
to black holes [30], and has furthermore been found for other classes of black hole solutions 
[3"H I3"2l [33"1 [3^] . In particular, this has the consequence that if we would live on M. A x S 1 
then from a four-dimensional point of view one would have an infinite non-uniqueness for 
static black holes of size similar to the size of the extra dimension, thus severely breaking the 
uniqueness of the Schwarzschild black hole. 

Another consequence of the new multi-black hole configurations of this paper is for the 
connection to uniform and non-uniforms strings on the cylinder. As mentioned above, there is 
evidence that the black hole on the cylinder phase merges with the non-uniform black string 
phase in a topology changing transition point. It follows from this that the copies of black 
hole on the cylinder solution merge with the copies of non-uniform black strings. However, 
with our new solutions, we add a continuous span of solutions connected to the copies of the 
black hole on the cylinder. Therefore, it is natural to ask whether the new solutions also 
merge with non-uniform black string solutions in a topology changing transition point. If 
so, it probes the question whether there exist, in addition to having new black hole on the 
cylinder solutions, also new non-uniform black string solutions. Thus, the new solutions of 
this paper presents a challenge for the current understanding of the phase diagram for black 
holes and strings on the cylinder. 

Another connection between strings and black holes on the cylinder is that a Gregory- 
Laflamme unstable uniform black string is believed to decay to a black hole on the cylinder 
(when the number of dimensions is less than the critical one p3]). However, the new solutions 
of this paper means that one can imagine them as intermediate steps in the decay. 

The solutions presented in this paper are clearly in an unstable equilibrium. Any small 
change in the position of one of the black holes on the cylinder will mean that the black holes 
will go even further out of balance, and the endpoint of this instability will presumably be a 
single black hole on the cylinder. Nevertheless, one can argue for their existence for example 
by imagining two equal size black holes on the cylinder, and then having mass thrown towards 
only one of the black holes in the same way from both sides of the black hole, i.e. that the 
solutions keep the inversion symmetry around both of the black holes. Then the matter will 
increase the size of one of the black holes, leaving the other of the same size. 

The construction of multi-black hole solutions also enables us to examine the possibility of 
further new types of black hole solutions in Kaluza-Klein spacetimes. In particular, analysis 
of the three-black hole configuration suggests the possibility that new static configurations 
may exist that consist of a lumpy black hole (i.e. 'peanut-like' shaped black objects), where 
the non-uniformities are supported by the gravitational stresses imposed by an external field. 
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The outline of this paper is as follows. In Section [2] we construct the new multi-black 
hole configurations on the cylinder to first order in the total mass of the system. In Section 
[3] the equilibrium condition for these configurations is explored, and a copying mechanism is 
presented that generates new equilibrium configurations from known ones. The first-order 
corrected thermodynamics of the multi-black hole solutions is given and analyzed in Section 
[31 We then present in Section \5\ the multi-black hole configurations in the phase diagram 
for Kaluza-Klein black holes, together with the already known black hole and black string 
solutions. Section [6J contains a more detailed analysis of the two simplest multi-black hole 
configurations, namely with two and three black holes. Finally, Section [7] contains a summary 
of our results, a discussion on its implications for possible new black hole and string phases 
and open problems. This concluding section also discusses in the context of an analogue 
fluid model a possible, but more speculative, relation of the multi-black hole configurations to 
configurations observed in the time evolution of fluid cylinders. Appendix A contains formulae 
that are used to compute thermodynamic quantities for the case of two unequal mass black 
holes on a cylinder. 



2 Construction of multi-black hole configurations on the cylin- 
der 

In this section we construct explicitly new solutions describing multi-black hole configurations 
on the cylinder, in the limit when the total mass of the black holes is small. 

2.1 General idea and starting point 

In the following we shall construct new solutions for multi-black hole configurations on the 
(f-dimensional cylinder x S 1 . The solutions are static and they describe configurations 

with several small black holes located at different points of the cylinder x S 1 . 

We require that all of the black holes are placed in the same point of the part of the 
cylinder. This is necessary in order to have equilibrium. Since all the black holes are placed 
in the same point of R rf_1 we can require the solution to be spherically symmetric on 
Since the solutions should solve the vacuum Einstein equations, the spherical symmetry has 
the consequence that we can write the metric for the multi-black hole configuration using the 
ansatz [TBI 1251135] 

ds 2 = -fdt 2 + jdR 2 + j^dv 2 + KR 2 dn 2 _ 2 , / = 1 - ||I , (2.1) 

where A(R, v) and K(R, v) are functions of the two coordinates R and v. As we shall see 
more explicitly below, the event horizons for the black holes are all placed at R = Rq. For 
simplicity, we set the radius of the cylinder to be 1. Thus, the R and v coordinates can 
be thought of as being measured in units of the radius of the cylinder. The v coordinate is 
periodic with period 2tt [16J. For R ^> 1, we are in the asymptotic region where the metric 
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asymptotes to the flat cylinder metric 

ds 2 = -dt 2 + dr 2 + r 2 dn 2 d _ 2 + dz 2 , (2.2) 

where z is periodic with period 2tt. Thus, we require that A(R, v ) — > 1 and K(R, v) — ► 1 for 
i? — > oo, and we see that R/r — > 1 and u/z — > 1 for it! — > oo. 

We construct in the following the metric for multi-black hole configurations on the cylinder 
M. d ~ l x S 1 in the limit where each of the black holes are small relatively to the distance between 
them. To this end, we employ the methods of [T7] to find the solution to leading order in the 
limit of small total mass. One can equivalently use the methods of \18\ [T9] to construct the 
metric. 

We proceed in the following to construct the solution in three steps: 

• Step 1: We find a metric corresponding to the Newtonian gravitational potential sourced 
by a configuration of small black holes on the cylinder. This metric is valid in the region 
R> R - 

• Step 2: We consider the Newtonian solution close to the sources, i.e. in the overlap 
region i?o<i?<l. 

• Step 3: We find a general solution near a given event horizon and match this solution 
to the metric in the overlap region found in Step 2. The resulting solution is valid in 
the region Rq < R <^ 1. 

With all these three steps implemented, we have a complete solution for all of the spacetime 
outside the event horizon. 

Note that the solutions that we find below generalize the previously studied case of a single 
black hole on a <i-dimensional cylinder |16 | I17 | HH 1 I19|. i.e. a black hole with S d ~ l topology in 
a d + 1 dimensional Kaluza-Klein space- time M d x S 1 , M d being d-dimensional Minkowski 
space. The solutions furthermore generalize the so-called copies of the single-black hole on 
the cylinder solution, corresponding to copying the solution several times across the cylinder, 
thus giving a multi-black hole solution where each of the black holes have the same mass and 
with the black holes placed equidistantly along the circle direction of the cylinder |28|, [29] . 

2.2 Step 1: The Newtonian region 

We construct here the linearized solution for the multi-black hole configuration in the region 
R^> Ro away from the event horizons. We require the black holes to be small such that they 
interact through Special Relativistic gravity (i.e. a Lorentz-invariant extension of Newtonian 
gravity). In such a Special Relativistic gravity theory we have a potential for each component 
of the energy-momentum tensor that we turn on. For static solutions on the cylinder it is 
well-known that the two relevant components of the energy-momentum tensor are the mass 
density g = Too and the binding energy (tension) b = —T zz [12]. These components source 
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the two gravitational potentials 

V 2 d> = 8.G N ^, V 2 5 = -^6, (2.3) 

where Gn is the (d + l)-dimensional Newton constant. From the components of the energy- 
momentum tensor one finds the total mass M and the relative binding energy (also known as 
the relative tension) n as |12j 

J d d x q(x) , n = — j d d x b(x) . (2.4) 



M 



In the limit of small total mass, we have that the relative binding energy goes to zero for a 
single black hole, i.e. n — ► for M — ► [T7]. From this we have that B/(G^M) — ► for 
M — > 0. Since is proportional to G^M, this means that we can neglect the binding energy 
potential B as compared to the mass density potential <£, since B goes like (G^M) 2 for small 
masses. With this, we see that we only need to consider the potential and we thus see that 
we are considering Newtonian gravity, with the only potential being the potential sourced 
by the mass density. 

We now proceed to find the Newtonian gravity potential <£. We consider a configuration 
of k black holes placed on the cylinder. We write M as the total mass of all of the black holes. 
Define v\ as the fraction of mass of the i th black hole, i.e. 

k 

M i = v i M, Yl v i = 1 > ( 2 - 5 ) 

where Mi is the mass of the i th black hole. Note that < v% < 1. As discussed above, we 
place the black holes in the same point of the B^ -1 part of the cylinder. This corresponds to 
r = in the (r, z) coordinates of the cylinder (|2.2|) . Let now z* be the z coordinate for the 
i th black hole with mass i/jM. We can then solve the equation for $ in (|2.3|) as 

*M) = - J^" M F(r,z), (2.6) 

with 

k oo 

F(r,z)=Y, E " SET" ( 2J ) 

i=l m=-oo V 2 + (z- z* - 2vrm) 2 ] 2 

The potential A2.6H thus describes the Newtonian gravitational potential sourced by our multi- 
black hole configuration. One can also write the function F(r, z) as the Fourier series 

kri ( k °° \ 

F ( r ' *) = ( 1 + 2 E v < E /l ( mr ) cos [ m ( z " ^*)] ) • (2-8) 

\ i=l m=l ' 



Here the constant kd is defined as 



1 d-2tt d _ 1 , . 
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and h(x) as 



d-5 1 d-3 

h(x) = 2 2 — X 2 K d _ 3 



(2.10) 



where /i(0) = 1, and K s (x) is the modified Bessel function of the second kind (in standard 
notation 1351). For r — > oo we see that 



F(r, z) 



r d—3 



(2.11) 



Inserting this in (12. 6p we verify that the potential <& has the correct asymptotic behavior for 
r — ► oo of a Newtonian potential on the cylinder describing an object with total mass M. 

We now proceed to find a metric in the form of the ansatz (|2.ip describing the linearized 
solution of the Einstein equations corresponding to the potential (|2.6|) . We first notice that in 
the ansatz (|2.ip we have that gu = — 1 + Rq 3 / R d ~ 3 ■ However, to leading order in G^M we 
have that gu = —1 — 2$. Therefore, we get that R~ d+3 is proportional to $(r, z). Demanding 
furthermore that Rjr for r —* oo, we are lead to define R as function of r and z as |16| 



R(r, z) 



F(r, z) 



i 

d-3 



(2.12) 



Thus, we see that in order for the linearized metric to fit into the ansatz (|2.ip . we need to 
define R as ()2. 12[) for the flat space metric. The choice of R ()2.12p is consistent with having 
the horizon at R = Rq since we see that defining R in terms of F(r, z) means that we are 
defining R to be constant on the equipotential surfaces of $ [16]. Since (]2.12|) defines R for 
the flat space metric, we need also to find a corresponding v(r,z) for the flat space limit of 
the ansatz (|2.ip . One can check, using the flat space metric (|2.2p in cylinder coordinates r 
and z, that in order to obtain a diagonal metric in the R and v coordinates, we need v to 
obey the partial differential equations |16| 



„<2-2 



„d-2 



Or 



-d z F(r,z) , 



0, 



-d r F(r, z) . 



(2.13) 



(d - 3)k d z " ' n (d- 3)k d 

Using the Fourier expansion (|2.8p of F(r, z) we find the following explicit solution for v(r, z) 

k +oo 



171 



-h(mr) 



d-3 



rh'(mr) 



(2.14) 



V = Z + 2 ^2 v i ^2 Sm [ m ( z — z i 
i=l m=l 

where h'(x) = dh(x)/dx. We see that u/z -> 1 as required above. Given the two coordinates 
R and v defined in (|2.12p and (|2.14p in terms of r and z, we can now find the corresponding 
flat space metric that can be written in the ansatz (|2.ip . We find the flat space metric 

An 



ds z 



-dt 2 + A dR 2 + 



K, 



J^dv 2 + K R 2 dQ 



d-2 ; 



(2.15) 



with the function Ao(r, z) and Ko(r,v) given by 

n d-2 



K (r,z) 



„2jL d-3 



F(r,z) 



(2.16) 
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Using now (|2.16p together with (|2.12p and (|2.14p . we can find the two functions Aq(R,v) and 
Kq[R, v) and we have thereby specified completely the flat space metric (I2.15p . 

With the flat space metric (|2. 15fl . as found above from requiring gu in the ansatz (|2,ip 
to be consistent with the Newtonian potential (|2,6p . we are now ready to find the complete 
metric to first order in G^M in the Newtonian regime R S> Rq. This problem is solved in 
general in [T7], and we refer to section 4.1 in that paper for the details. The upshot is that 
given the flat space metric (|2,15p defined from the Newtonian potential $ in (|2.6p . we can 
find the correction to first order in Gt$M of the functions A(R, v) and K (r, v) as 

i 1 1 - , • , ,„•;::_. i t„ - ,,. r „^_\ a , 




Thus, given Aq(R,v) and Kq(R,v), as found above in (|2.16p . <|2.12j> and <|2.14j> . we can find 
A(R, v) and K(R, v) to first order in G^M, or, equivalently, to first order in -Rq -3 - Combining 
this with the ansatz for the metric (12. IB . we have actually found the metric up to first order 



in Rq 3 (i.e. in M) in the Newtonian region R S> -Rq, for any given distribution of k small 
black holes on the cylinder. 

2.3 Step 2: The overlap region 

In the previous Section 12.21 we found the metric for any given distribution of k small black 
holes on the cylinder to first order in the total mass. This metric is valid for R 3> Rq, i.e. away 
from the horizon. In this section we examine now this solution in the region Rq <C R <C 1, 
which we dub the overlap region, since this is the region where both the Newtonian regime 
and the near- horizon solutions are valid. As we shall see below, the analysis of the solution 
in the overlap region gives in turn a restriction on what configurations of black holes that we 
can find a metric for, namely that the k black holes should be in equilibrium with each other 
with respect to the Newtonian gravitational forces between them. 

Before turning to the first-order corrected metric found in Section 12.21 we first consider 
how the potential $ looks when going near the sources, and subsequently how the flat space 
metric (I2.15P behaves. In terms of the flat space coordinates R and v found in (|2.12p and 
(|2.14p . this corresponds to having R <C 1. Note that since we have k small black holes we 
have to specify to which of these we are close. In line with this, it is useful to define for the 
i th black hole the spherical coordinates p and 9 by 

r = psinO , z — z* = p cos 6. (2-18) 

Notice here that the angle 6 is defined in the interval [0, 7r]. We then conclude from (|2.18p 
that going near the i th black hole corresponds to having p 1. We begin by examining the 
function F(r, z) in (|2.T|) near the i th black hole. In terms of the spherical coordinates (|2.18p 
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we find that 



F(p, 9) = Vip-^-V + A« + Af } cos 6 p + A^(dcos 2 - l)/? 2 + O (p 3 ) , (2.19) 
for p <C 1, where 



A «) 



2((d-2) 

W-2 



+E 



(2tt)- - — I ■• 



d-2 



(27T) 



d— 2 



C ( d - 2 ' 1 -^ +C l d - 2 ' 1+ 27 r 



,(2.20) 



A«=(d-2)x; 



i=i 



-1 



(2vr 



id— 1 



(2vr) c 



Here 



C(«-M-2)-C(rf-i,l + £ 



2 ^ (2tt) 
i=i L^j v ; 



C(,,i-|) +C ( rf ,i + | 



(2.21) 



(2.22) 



C(s,l + a) = + 

m=l 



m + a ^ , 



(2.23) 



is the Generalized Riemann Zeta function and labels the distance in the z direction between 
nth „_ j „-th 



the j and i black hole as follows 



if < z* - z* < 2vr , 



Zij = 2ir + Zj 



if 



2vr < z* - z* < . 



(2.24) 



We see that this definition ensures that < < 2ir. 



Using now (|2.19p - (|2.24p with (|2.6p one obtains the behavior of the Newtonian potential 
near the z th black hole. 

From the potential for p <C 1 obtained by inserting (|2.19|) in (|2.6p we see that the 
first term in (|2.19p corresponds to the flat space gravitational potential due to the i th mass 
Mi = v%M and the second term is a constant potential due to its images and the presence of 
the other masses and their images! 2 ] Furthermore the third term in (|2.19p is proportional to 
/ocos 6 = z — z* and therefore this term gives a non-zero constant term in <9 Z <3? if we have that 
A j 1 given in (|2.21[) is non-zero. This therefore corresponds to the external force on the i th 
black hole, due to the other k — 1 black holes. In Section [3] we verify this interpretation. 

Since A^ is proportional to the external force on the i th black hole, it is clear that one 



cannot expect a static solution to exists if A^ is non-zero, since then the i black hole would 

a (i) 

accelerate along the z axis. Therefore, the only hope of getting a static solution is if A\ =0 
for all i = 1,2, k, i.e. that the external forces on each of the k black holes are zero. When 



2 In particular, the origin of the three terms contributing to A' 1 ' in (|2.20[) is as follows. The first term comes 
from the images of the i black hole, the second term from the other k — 1 black holes and the third term 
from the images of these. 
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constructing our solution, we therefore assume that = for all i. From (|2~2T|) . wc sec 
that this gives conditions on the relation between the positions z* and the mass ratios z/j. 
We explore these conditions further in Section [3J Note that the equilibrium established with 
kf = for all i is an unstable equilibrium, i.e. generic small disturbances in the position of 
one of the black holes will disturb the balance of the configuration and result in the merger 
of all of the black holes into a single black hole. 

We consider now how the flat space metric fj2. 15f> looks near the black holes. To this end, 
it is useful to consider the flat space coordinates R and v found in (|2,12p and (|2.14p near the 
i th black hole. Using (|2,19p . we see that 

R d - 3 ^vT 1 k d p d - 2 , v ~ Vl - Vl ^fl k f f dx(smx) d - 2 , (2.25) 

a — 6 J x=0 

for p <C 1, with the number pi defined as 

i-l 

Pi = it , for i = 1 ; pi = 7T — 2n Vj , for i = 2, • • • , k . (2.26) 

3=1 

Note that 9 = corresponds to v = pi and = ir corresponds to v = pi — lixvi. So the range 
of the coordinate v can belong to one of the k intervals li defined as 

k 

h = [Pi - 2nvi,pi] , with (J li = [-7T, 7r] , (2.27) 

i=l 

where the last condition follows from the fact that Yli=i v % = ^-- The physical meaning of the 
intervals (|2.27[) is that each of the interval corresponds to one of the black holes. So, being 
close to the i th black hole in (R, v) coordinates corresponds to having R <^ 1 and v G li. This 
feature continues to hold also in the first-order corrected metric. 

In order to match the metric in the overlap region to the metric near the horizons of 
the black holes, it is natural to change the ansatz (|2,ip into a form which resembles more 
the spherical coordinates (p,0), instead of the cylindrical coordinates (r, z). Given a solution 
in the form of the ansatz (|2.1|) with the functions A(R,v) and K(R,v), we define therefore, 

hole, the new coordinates p and 

— 1 i ~d— i d 2 , _ 

a k d p , v=pi-Ui- -k d 

a — o J x=0 

where pi is defined in (|2.26p . and 9 = corresponds to v = pi while 8 = tt corresponds to 
v = Pi — 2nui. The coordinates (p,6) are defined such that p = p{R) and = 9{v) and such 
that for the flat space metric we have p ~ p and 6 ~ 9 for p <C 1, as one can see from ()2.25p . 
We define furthermore the two functions A(p, 9) and K{p, 9) by 

A = ^^(u~ 1 k d p)~^ A, K = sin 2 9{u^ k d p)~^ K, (2.29) 



d K(i 



relative to the i th black hole, the new coordinates p and 9 by 

R d - 3 = „r 1 k d p d - 2 , v = pi - Ui^f^k, 1 I dx{smx) d - 2 . (2.28) 



3 The factor Vi in the second expression of (|2.28|l guarantees that 6 — <-> v — Pi while 9 — n <-> v = Pi — 2iWi 
We choose to include the Vi in the first line of (|2.28jl to have p/p — > 1 when R — > 0; see (12.3411 . 
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and the parameter pr> by 

JL-2 _ k -l r>d-3 
Po — K d a ' 

such that we can write the ansatz (|2.ip in the alternative form 

A 



ds l 



-fdt 2 + jdp 2 + -^- 2 P 2 d6 2 + Kp 2 sin 2 0dfi 2 _ 2 , / = 1 



p^ 2 



(2.30) 



(2.31) 



Note that the event horizon for the i th black hole is located at p = uf~ 2 p®. 

Turning to the flat space metric, corresponding to the zero total mass limit of the metric 
for the multi-black hole configuration, we can reformulate the above results for the (R, v) 
coordinates in terms of the (p, 9) coordinates. We write the flat space limit of the ansatz 



ds 2 = -dt 2 + A dp 2 + P 2 ^ 2 + Kop 2 sin 2 9dn 2 d _ 2 



(2.32) 



The functions Aq(p, 9) and Ko(p, 9) defining the flat space metric (|2.32j) are most easily found 
using the relations 



Ac 



(d p p) 2 + p 2 K 



p 2 sin 2 i 
p 2 sin 2 i 



(2.33) 



Implementing now the definitions (|2.28p and the results (|2.25j) . we see that for p <C 1 (which 
is equivalent to p <C 1) we get the expansion^ 



P = P 



sm 



1 + 



sm 



d-2 



pd-2 + ^Jh_ cos6p d-i 



i + 



d-2 
ia(0 



+ o(~ P d ] 

cosf/r - + 0(j> 



(d- l)(d-2) 

Using this with (12.331) . we find the following expansions for Ar,(p,9) and Ko(p,0) 



An 



1 + 



^^-W), „_ 2 



1 + 



+ C(p d 



(2.34) 



(2.35) 



for p <C 1. We included here the corrections up to order p d 2 . Note that the next corrections 
come in at order p d since here and in the following we have set A^ = 0. 

Having understood the flat space metric in the ansatz (|2.3ip near the i th black hole, we 
are now ready to collect all the results and write down a first-order corrected metric near the 
i th black hole. First, we note that using the definition (|2.29p it follows from the general form 
(|2.17|) for the first-order corrected metric in the (R, v) coordinates that we obtain the general 
form for the first-order corrected metric in the (p, 9) coordinates, 



A = Ac 



P 



^pV 



2(d - 2) p d - 2 & ' pA ° ' 



K = K n 



P 



^pV 



2(d - 2) pd-^ Ko ■ 



(2.36) 



4 We included here for completeness the A[ terms although we set A]" = in the actual solutions in order 
to have a static solution, as discussed above. 



(i) 
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Given the full flat space functions Aq(p,9) and Ko(p,9), the functions A(p,9) and K(p,9) 

in (I2.36|) when inserted in the ansatz (I2.3ip describe the first-order corrected metric for a 

configuration of small black holes in the region p S> pq. Using now the p C 1 expansion of 

Aq an d Kq found in (|2.35p we get the following explicit expansions of the first-order corrected 
i 

metric for v^ 2 p <C p <C 1 



A ~ 1 + 



Tp^-v lP f , #-! + -!__ 1~p a - l -v lP t l . (2.37) 



d-2 

Thus, the functions (|2,37p with the ansatz (|2.3ip give the metric of the multi-black hole 

i 

configuration in the overlap region v^~ 2 po <p<l. In Section \2A\ we shall match this with 
the metric in the near-horizon region. 

Regularity of the solution 

We can now address the regularity of the multi-black hole solution given the above results for 
the first order correction. We already argued above that we need the equilibrium condition 
A? = to hold for all i = 1, k, since otherwise the configuration that we are describing 
cannot be static. However, this should also follow from demanding regularity of the solution, 
since with a non-zero Newtonian force present on the black hole the only way to keep it static 
is to introduce a counter-balancing force supported by a singularity. Therefore, it is important 
to examine the regularity of the solution corresponding to (|2,36p with or without the presence 
of the terms. 

For a metric in the form of the ansatz (|2.3ip . one can have singularities for 9 — > 0, tt, since 
the metric component along the (d — 2)-sphere goes to zero there. A necessary condition to 
avoid such singularities is that for 9 — > 0, tt the 9 part plus the (d — 2)-sphere part of the 
metric (|2.3ip becomes locally like the metric of a (d — l)-sphere d9 2 + sin 2 9dQ^_ 2 since then 
9 = 0, tt corresponds to the poles of the {d — l)-sphere. This is only the case provided that 

A 

1 for 9 -> 0, tt . (2.38) 



K d-i 



Therefore, we should examine under which conditions the correction (|2.36p obeys Eq. (|2.38p . 
First, let us assume that the flat space functions Aq, Kq obey Eq. (|2.38p . «-e. Aq/Kq -1 — > 1 
for 9 — > 0, 7T. From this one can infer that dp log Aq — (d — l)dp log Kq — > for 9 — > 0, tt. Using 
this, it is not hard to check that Eq. (|2.38p is fulfilled with A and K given by (|2.36p . Thus, in 
order to fulfil (|2,38p we only need to check that it is fulfilled for the flat space metric. This is 
indeed found to be the case, both for the A® terms and the A^ terms. Thus, the metric is 
regular at the poles 9 = 0,tt also with the external force on the z th black hole being present. 
This is presumably because we cannot see the irregularity of the solution at this order since we 
can neglect the binding energy, which accounts for the self-interaction of the solution. Thus, 
we expect singularities to appear at second order in the total mass for solutions which do not 
obey the equilibrium condition Aj = 0. 
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2.4 Step 3: The near-horizon region 

In Section 12.21 we found the metric (to first order in the mass) for a general multi-black hole 
configuration in the Newtonian region R ^> Rq. We now complete the metric for the multi- 
black hole configuration by finding the metric near the horizon. This is done by matching 
with the metric in the overlap region Rq -C R <C 1, as found in Section [2.31 

Take the metric (|2.3ip with (|2.37p which describes the geometry near the i th black hole, 

i 

i.e. in the overlap region vf~ 2 po <C p <C 1. We notice here the key point that A and K are 

independent of 9. This means that we can assume that A and K are independent of 6 for 

i 

v^~ 2 Po <p<l. The next step is therefore to find the most general solution of the vacuum 
Einstein equations for a metric of the form (|2.31|) with A = A{p) and K = K{p), i.e. without 
any 6 dependence. This gives the result |17j 

d-2 _ jo 1 — 1/J 2 2 

A'^J =K~^ = W p +w , (2.39) 

w ViPl 

where w is an arbitrary constant. Note that, setting w = 1, the ansatz (|2.31j) with (|2.39p 
describes the (d + l)-dimensional Schwarzschild black hole solution. 

We can now fix this constant w by matching the functions (|2.39p to the behavior of A and 
K in the overlap region (|2.37|) . This yields 

-l + yP^ + ^n. (2-40) 

Thus, using (|2.39p with (|2.40p in the ansatz (|2.31j) . we have obtained the metric for a general 

multi-black hole configuration, in the limit of small total mass, in the near-horizon region 
i 

vf~ 2 po <p<l. Supplementing this with the metric in the Newtonian region R^> Rq found 
in Section I2.2[ we see that we have obtained the full metric for the general multi-black hole 
configuration to first order in the mass in the limit of small total mass. 

Inserting (|2.39p and (|2.40p in the ansatz (|2.31j) . we can write the near-horizon metric near 
the i th black holes located at (r, z) = (0, z*) as 

ds 2 = -fdt 2 + /- 1 G" 2 ^ i dp 2 + G~^p 2 (d6 2 + sin 2 9 d^_ 2 ) , (2.41) 

where (up to first order in Pq~ 2 ) 



P uj Uip Q 4 

1 

The horizon is located at p = vf~ 2 po and the range of 9 is from to n. 

3 Equilibrium configurations 

From the results of Section [2] we have that near the i th black hole the gradient of the gravi- 
tational potential along the z-direction is 

= 7-^- hi - - Af» + owl , (3.1) 



(d-l)fl 



d-l 
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for p <C 1. The first term is evidently the gravitational attraction due to the mass of the i th 
black hole, while the second term is a net force on the i th black hole, which originates from 
the other k — 1 black holes and their images in the configuration^ Having such a force on 
the i th black hole is clearly not consistent with having a static solution. Therefore, as already 
discussed in Section 12.31 we require that the solutions fulfil the equilibrium condition 

Aj =0 for i = 1, (3.2) 

In Section 13.11 we explore this condition further, and we describe a method of how to find 
configurations, i.e. a set of masses V{ and positions z*, such that the equilibrium condition 
(|3.2|) is fulfilled. We furthermore describe in Section 13.21 how to generate new equilibrium 
configurations from known ones by copying. 

As already discussed in Section [2T3"1 the equilibrium of the k black holes is unstable towards 
perturbations in the positions of the black holes. We compare this physical intuition with the 
results for the two-black hole solution in Section T6. 11 

3.1 Construction of equilibrium configurations 

In the following we describe a construction method that allows one to find equilibrium con- 
figurations fulfilling (|3.2p . While doing so we further clarify the equilibrium conditions. 

Condition f|3.2|) per se is not in general sufficient to identify specific parameters of config- 
urations that are in equilibrium. In the following we describe a procedure from which we can 
obtain an equilibrium configuration given a set of black hole positions (with some restrictions). 

(i) 

We first note that we can write A\ as a sum of the potential gradients corresponding to 
the gravitational force due to each of the k — 1 other black holes on the i th black hole aqj 

k 

where Vij corresponds to the gravitational field on the i th black hole from the j th black hole, 
given by 



^ = (rf - 2 )(J IT -^_L_[ C ( d - M _g)_ ( ( d _ 1 , 1 + i ) >, (3 .4) 



for j ^ i. We can now furthermore define Fij = ViVjVij as the Newtonian force on the i th mass 
due to the j th mass (and its images as seen in the covering space of the circle). Of course, 
to obtain the actual Newtonian force we have to multiply F{j with 8ttG^M 2 /{{d — l)Qd-i)- 
With this, we can write (|3.3p as the condition of zero external force on each of the k masses 

k 

£ Fij=0, (3.5) 



5 The images of the i th black hole only contribute in Eq. (|3.1[l in the terms of O(p). 
6 Note that the force on the i th black hole is A^' 8ttGn M / ((d - 
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for % = 1, k. We can now verify an important property, namely that Newton's law F\j — 
—Fji is satisfied. Clearly this is equivalent to Vij = —Vji. From (13.41) and the definition (|2.24l) 
of Zij for the i th black hole, we see that Vij = —Vji follows from the following identify for the 
Generalized Zeta function (|2.23p 

fey^(^) + <(^-i)-(T) S + c(M-a-c(M + i).(3,) 

We now illustrate our procedure of finding equilibrium configurations by considering the 
k = 3 black hole case. The generalization to an arbitrary number of black holes is easily done. 
First, consider a given set of positions of the black holes (z-[,z 2 ,z 3 ). From these positions 
we get from (|3,4p . We now want to find u\, u 2 and u 3 such that we get an equilibrium 
configuration. From f)3.5j) we see using Fij = —Fji that there are only two independent 
equations, which we can write as v 2 Vi 2 + v 3 V\ 3 = and —V\V\ 2 + ^3^23 = 0. Using now that 
v 3 = 1 — v\ — v 2 , we get the following result for ui, v 2 and v 3 

_ ^23 Vl3 _ V[2 , „^ 

Vl ~ v 12 - v l3 + v 23 ' 1/2 " v 12 - v l3 + v 23 ' 1/3 " v l2 - v 13 + v 23 ■ [ ■ ) 

Thus, we see that choosing the positions of the three black holes gives us Vij which again 
gives us z^i and v 2 from (|3.7p . implementing the zero force condition (|3.5p . 

However, it is important to note that we need to impose the physical requirement of having 
only positive masses, i.e. < < 1 for all i. This again gives restrictions on the positions 
that one can choose. For k = 3 one can check that these restrictions are satisfied under the 
fairly mild conditions z\ = < z 2 < tt < z\ < 2ir and z^ — z 2 < w. 

The above construction method that we described for k = 3 can be extended to con- 
figurations with any number of black holes subjected to some constraints on their relative 
positions. One then solves the k — 1 independent zero force conditions from (|3.5p for the 
k — 1 independent mass parameters z/j. Note that one can infer from this way of solving the 
equilibrium condition (|3.2|) that in general a k black hole configuration has k independent 
parameters, e.g. the rescaled mass and the k — 1 positions^ Another way to see that we have 
k independent parameters for a configuration with k black holes is to note that by specifying 
that z* < z 2 < ... < z], and by giving the k absolute masses z^M (or alternatively the rescaled 
total mass and k — 1 of the mass parameters Vi) we can determine an equilibrium configuration 
using the analysis above. 

3.2 New equilibrium configurations by copying 

We described above a general method to build equilibrium configurations. In this section we 
consider a way to generate new equilibrium configurations using already known ones. This is 

7 Note that there are special configurations with a high amount of symmetry where the mass ratios Vi are not 
fixed given the positions z* . An example of this is the two black hole case with z\ = and zj = 7r. However, 
the number of independent parameters is always k for a k black hole configuration, i.e. for the two black hole 
case the two parameters can be taken to be fj, and v\. 



15 



done by copying the configurations a number of times around the circle. This generalizes the 
copies of the single-black hole solution (28J EE E] . 

We imagine a configuration given with k black holes, specified with the positions z* and 
masses Ui, i = 1, ...,k. We assume this configuration is in equilibrium, i.e. that (|3,2p is 
satisfied. We also assume that the positions are ordered such that < z* < z* +1 < 2-it for 
i = 1, k — 1. Given now an integer q, we can copy this configuration q times, to obtain a 
new equilibrium configuration as follows. We define^ 

z* +nk = -(z* + 2vrn) , u i+nk = -z/j , (3.8) 

for % = 1, ...,k and n = 0, ...,q — 1. Then zi, ■ ■■,z kq and ...,Vkq defines a new configuration 
with kq black holes. In particular we have that Yla=i = 1 an d that < z* < z* +1 < 2-k for 
a = 1, kq — 1. 

We first verify that the new configuration of kq black holes obeys the equilibrium conditions 
(]3.2p . Note that this check is needed only for the first k black holes (out of the kq black 
holes) since the black hole configuration is symmetric under the transformation z* — > z* +k , 
v a ^a+k if we furthermore make the displacement z — > z + 2ir/q. Consider therefore the 
zero force condition on the i th black hole, with i = 1, k. Using f|3.3|) we can write this as 



q-l q-l k 

Vi+nkYi,i+nk + ^j+nkVij+nk = , (3-9) 

n=l n=0j=l,jj=i 

with V a b given by (|3.4|) . Here we have split up the contributions such that the first term 
corresponds to the copies of the i th black hole, while the second term corresponds to the other 
k — 1 black holes and their copies. Using now that Zij+nk = ~p an d ^i,j+nk = -f- + 
as one can infer from the definition (|2.24p . it is straightforward to verify, with the aid of the 
definition of the generalized Zeta function (|2.23|) , that we have 

q-l q-l 

V i)i+nk = , J2 ViJ+nk = q d ~% ■ (3.10) 

n=l n=0 

Using this, we see that it follows from the equilibrium condition Y2j=i j^i v jVij = f° r the k 
black hole configuration that the equilibrium condition (|3.9p is satisfied for the kq black hole 
configuration. 

It is useful to consider how one can express the metric for the q copied configuration in 
terms of the metric for the k black hole configuration. To this end, we note that one easily 
sees from (|2.7p that 

F(r, z) = q d - 3 F(qr, qz) . (3.11) 

This gives in turn that Ao(r, z) = Ao(qr,qz) and Ko(r,z) = K${qr,qz). By carefully using 
these relations, we infer that Aq(R,v) = Ao(qR,qv) and Kq(R,v) = Ko(qR,qv). Therefore, 



Note that here and in the following we put a hat symbol on all the functions, parameters and quantities 
that correspond to the new configuration that we copied q times. 
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we have from (|2.17p that 



A(R,v)=A(qR,qv), K(R,v) = K(qR,qv) . (3.12) 

From this we can read off the metric for the q copied configuration in terms of the metric 
for the k black hole configuration. Notice that this relation precisely corresponds to the one 
found in |29j from a more general point of view. 



4 Thermodynamics of the multi-black hole configuration 

In this section we begin by determining the thermodynamic properties of the multi-black hole 
configurations. This is accomplished in Section [4.11 We subsequently find in Section [4.21 that 
the obtained thermodynamics is consistent with a simple Newtonian interpretation. 

4.1 Thermodynamic properties 

In this section we find the thermodynamic quantities for multi-black hole configurations on 
the cylinder to first order in the mass in the limit of small total mass. 

We begin by considering the quantities that one can read off from the event horizons. For 
the i th black hole the metric near the horizon is given by (|2.4ip - (|2.42p . The temperature is 
now found in the standard way by computing the surface gravity while the entropy is found 
from computing the area of the event horizon divided with 4Gn- This yields the following 
entropy Si and temperature Tj for the i th black hole 



with as defined in (|23U]h 

Turning to the asymptotic quantities, we need to determine the total mass M and the 
relative tension (binding energy) n. To determine M and n, we first notice the fact that the 
multi-black hole solution obeys the first law of thermodynamics [29| 



5M = ^TiSSi. (4.3) 



i=i 



This is derived in [29] using the ansatz (12. lj) for a single connected horizon, but the argument 
there is easily generalized to k connected horizons. Note that in (|4,3p we do not have the 
variation of the circumference of the cylinder since we have fixed the circumference to be 2n. 
This term is however easily added (see Ref. and below). 
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It is a general property of the ansatz (|2.ip that Rq, M and n are related as [29 



M= n d -2 pd - 3 (rf-i)(rf-3) 

8G N d - 2 - n 



(4.4) 



This is easily seen from considering the metric (|2.ip for i? — > oo. Using the definition of /?o in 
(j2.30p we can write this as 



M 



(4.5) 



16vrG N u d-2-n 

We can now insert (|4.ip . (|4.2|) and (j4.5|) into the first law (14. 3p for a given variation of po, 
which yields the following result 

k 

(4.6) 



d-2 



d — 2 5po 



8=1 



to first order in p . We used here that n — > for pg ~ > 0- From (|4.5|) and (|4.6|) we then 



conclude that M and n, to first order in p , are 



(d— l)fL_1 rl 9 



2(d-2)> 



1=1 



(4.7) 



i=l 



2(<2-2) N 



(4.8) 



Thus, the physical quantities relevant for the thermodynamics of the k black hole configuration 
are given by (|4TT|) . K2\\ . K7\> and (|Oj> . 

We consider now how the relative tension n and the entropies Si, as given above, behave 
as a function of the total mass M. To this end, it is useful to define the rescaled mass p as 



16ttG n M 16ttG n M 



L d~2 (27r) d - 2 ' 

where we used that the circumference L = 2n. Using now (|4.7p and 
function of /i is given by 



4(d- 1)J2 



(4.9) 

we get that n as 
(4.10) 



i=i 



We use this expression in Section [5] since it gives the linear slope of the multi-black hole 
configuration in the (p, n) phase diagram. There we also provide a rough estimate of the 
range of p for which (|4.10p is a good approximation. 

Turning to the entropies, we have that the entropy of the i th black hole, in terms of the 
rescaled total mass p, is 

d-2 / -. k 



Si(p) 



{2n) d - l {u i p)^ 
4n£*(d-l)^G N 



1 + 



(27T) 



2(d-2)Q d _ 



8=1 



.(4.11) 
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One can now compute the total entropy S , to tai(^) as the sum of the entropies (|4.1ip for each 
of the k black holes. 

As already mentioned, the k black hole configurations are unstable with respect to small 
changes in the positions of the black holes. Generic disturbances will destabilize the config- 
uration and presumably the k black holes will merge into a single black hole. Therefore, we 
expect in general that the entropy for a single black hole is always greater than the total 
entropy of the k black holes, for same total mass p, i.e. S to ta,\(p) < S'ibh(m)- This can indeed 
be verified from Eq. (|4.1ip . for sufficiently small p. We examine these questions in detail in 
Section [6.11 for the two black hole case. 

It is important to note that from the temperatures (|4.2|) one can see that they in general 
are not equal for the black holes in the configuration. This means that generically the multi- 
black hole configurations are not in thermal equilibrium. In fact, it is easy to see from (j4.2|) 
that the only configurations at this order that are in thermal equilibrium are the copies of the 
single-black hole solution studied previously in [281 ESI E! • 

4.2 Newtonian interpretation of the thermodynamics 

The variable po was useful to construct the multi-black hole solution but is not the most 
appropriate one for the physical interpretation of the solution and its thermodynamic quanti- 
ties (|4.ip - (|4.8p . since it does not have an invariant meaning. A more natural variable for the 
physical interpretation, as will be confirmed below, is the "areal" radius. We define a set of 
k "areal" radii poa), i = 1, • • • > k, by 

^ = vf-* po h+-^—-^\ . (4.12) 

Using this definition the first-order corrected horizon area of the i th black hole takes the 
appropriate form 

A® = JVipg , (4.13) 

for a (d — l)-sphere of radius pom. We can now rewrite, to leading order, the corrected 
thermodynamic quantities (I4.ip - (l4.8p in terms of these "areal" radii. 

The corrected entropy (|4,ip and temperature (|4,2p of the i th black hole takes the form 

o "d-l 

Si = Ar m , T, = T 0W (1 + <!>,), T 0(4) = — — , (4.14) 
4G N 47rp 0(i) 

where we have defined the potential 

A® 

From the form of A® in Eq. (f2^20|) we see that <I>j is precisely the Newtonian potential created 
by all images of the z th black hole as well as all other k — 1 masses (and their images) as seen 
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from the location of the z th black hole. The interpretation of the form for the temperature 
in (|4.14p is that T u\ is the intrinsic temperature of the i th black hole, i.e. when it would be 
isolated in flat empty (d+ l)-dimensional space. The second term is the redshift contribution 
coming from the gravitational potential ^j. 

Similarly, the total mass (|4.7p of the configuration can be written to leading order as 



k 
i=l 



M °« = IGttGn ' (416) 



where is defined in ()4.15|) . Again, the physical interpretation can be clarified as follows: 
The first term 



M = 2^M 0(i) , (4.17) 

i=i 

is the some of the individual masses M u\ when they would be isolated, while the second term 

1 k 

^Newton = j E M °« $i ' ( 418 ) 
i=l 

is precisely the negative gravitational (Newtonian) potential energy that appears as a conse- 
quence of the black holes and their images. 

From the above results it follows that one can derive the formula for the relative tension 
in (|4.8p by a purely Newtonian argument, as was first done for the single black hole case in 
Ref. [19J. To see this, note that when we also allow for the length L of the circle to vary, the 
generalized first law of thermodynamics (|4.3p reads 

5M = J2 TiSSi + ^j-5L , (4.19) 

i=l 

since T = nM/L is the tension in the circle direction. The relative tension can thus be 
computed from 

L (dM\ 

(4 ' 20> 

As described above, from a purely Newtonian analysis one knows that the total mass M = 
Mq + J/Newton is the sum of the intrinsic mass plus the gravitational potential energy given in 
(|4.18j) . Furthermore, the condition of keeping Si fixed means that we should keep fixed the 
mass Mg(i) of each black hole, and hence also Mq. It thus follows from (|4.20p that to leading 
order 

"-Sftl-ar-J^— «5£^ £ -5r ■ (4 - 21> 

where we used ()4.18j) . ()4.15p in the second step. To compute the derivative we need to know 
how A^) scales with L. While the expression for AW in (UD is for our choice L = 2ir, it is 
not difficult to see that keeping L arbitrary amounts to the rescaling A® — > (2ir / L) d ~ 2 A^. 
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Using this in (|4.21|) along with Mqu^/Mq = Vi immediately shows that we recover our result 
(|4.8p for the relative tension. 

As a consequence, we conclude that the entire thermodynamics of the first-order corrected 
multi black-hole solutions can be appropriately interpreted from a Newtonian point of view. 



5 Phase diagram for the multi-black hole configurations 

As mentioned in the Introduction, the whole set of different multi-black hole configurations 
are part of a larger set of black holes, black strings and other black objects which are asymp- 
totically M d x S l [21 0]. For this reason, it is very useful to depict the multi-black hole 
configurations in a (/i, n) phase diagram [12, 29J in order to understand the phase structure 
of all the solutions asymptoting to Ai d x S 1 . 

A multi-black hole configuration corresponds to a point in the (//, n) phase diagram. The 
coordinates of this point are given by (|4.7p - (|4.9p . However, since we look at the limit of small 
gravitational interactions, it is useful to have n as function of fj,. This is given by (|4.10p . 
Therefore, n(/j,) as given in (|4.10p is valid for small [/,. For a fixed n, one can then consider 
the range of n for a configuration with k black holes. This can be seen using the following 
inequality for a k black hole configuration 

2 C(d-2) * (i) d _ 3 2C(d-2) 

(27r) d-2 ~ k (27^-2 ' ^ 

where is defined in (|2.20p . The lower bound corresponds here to the single-black hole 
case (k = 1) while the upper bound corresponds to the case of k black holes of equal mass, 
distributed equidistantly around the cylinder. These are the so-called copies of the single 
black hole on the cylinder considered in [281 ISU HZ]- Now, using the inequality (|5.ip with 
()4.10p we see that in the (//, n) phase diagram the k black hole configurations corresponds to 
the points lying above the single-black hole phase and below the k copied phase. We have 
depicted this for d = 5 in Fig. Q]in the case of two black holes on the cylinder {k = 2). We 
see that the phases with two unequal black holes lie in between the single localized black hole 
phase (LBH) and the phase with two equal size black holes (LBH2). We have depicted here 
the phases using the numerical data found in |24j for the single black hole phase (LBH). Note 
that it is not clear that the phases of the two black hole configurations will stay in between 
the LBH and LBH2 phases when we go beyond our perturbative solution (see discussion in 
Section [7.2|) . In Fig. [I] we have furthermore depicted the uniform black string phase (UBS), 
which has n = l/(d — 2), and the non- uniform black string phase (NUBS), along with the 
two-copied non- uniform black string phase (NUBS2)I§ 

From (|4,10p and the inequality (|5.ip we see that for a given mass \i we have a continu- 
ously infinite non-uniqueness of solutions with k black holes. However, the non-uniqueness 



9 For the d = 5 non-uniform black string we have used the data given in [11U12] , The map to the two-copied 
solution is given in |29| . 
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Figure 1: Phase diagram for d = 5 with n versus [x for the two-black hole configurations 
spanning the area in between the single black hole (LBH) and two equal size black holes 
(LBH2). Moreover, we have drawn the uniform black string phase (UBS), the non-uniform 
black string phase (NUBS) and its two-copied phase (NUBS2). 

of solutions is even worse than this. If we consider a k black hole solution it is described by 
k parameters, as explained in Section 13.11 Thus, since the solutions with k black holes span 
a two-dimensional area in the (//, n) diagram we need k — 2 extra parameters, beyond /i and 
n, to point to a specific solution with k black holes. Therefore, there is a continuous infinite 
non-uniqueness of solutions for certain points in the phase diagram, when k > 3. Moreover, if 
we do not specify k but instead consider all possible multi-black hole configurations, we have 
an infinite layer of solutions in the phase diagram, since one can always consider adding a 
small black hole to a given multi-black hole configuration. 

Hence our results show a continuous non-uniqueness for solutions with fixed M. Such 
non-uniqueness was also observed in Ref. [30] for bubble-black hole sequences, which are not 
spherically symmetric on and lie in the region < n < d — 2 of the (//, n) phase 

diagram. The multi-black hole configurations of this paper are therefore the first example of 
continuous non- uniqueness for solutions that are spherically symmetric on R^ -1 . 

Considering the phase diagram for the two black hole configurations depicted in Fig. [H 
it is interesting to consider what happens when moving up in n. One way to do this is to 
increase \x such that the ratios are fixed. In this case the two black holes are growing and 
eventually their horizons will meet. Thus, the natural question is then what happens when 
approaching this point. There seems to be two possibilities: 

1) When the horizons of the two black holes meet, their temperatures are not equal, and 
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the solution will be singular in the meeting point. 

2) The temperatures of the two black holes will approach each other and when the two 
black holes meet they will merge into a new non-uniform black string phase different 
from both the original non-uniform black string phase emanating from the Gregory- 
Laflamme point, and the two-copied non-uniform black string phase. 

We explore these scenarios further in Section 16.11 In Section [7] we discuss the possible im- 
plications for the Gregory-Laflamme instability if there should exist new non-uniform string 
phases. 

Finally, we note that it is useful to give a rough estimate of the validity of the perturbative 
k black hole solutions found in Section [2j For this purpose we can employ the estimate made 
for the single black hole solution in [17] ■ A lower estimate can be found by considering the k 
copied phase, since we expect this to be the first solution for which the first order correction 
becomes invalid, as one increases p. We therefore take the function F(p, 9) in (|2,19p and 
consider when the contribution from the term with is equal to the one with AW. This 
happpens for p 2 ~ 8ir 2 ((d — 2)/(k 2 (d — l)(d — 2)£(d)). This can be used to get an upper bound 
for the Schwarzschild radius A; _1 /( d_2 ) p . Plugging that into p in terms of po, one obtains 
a rough upper bound on p. For k black holes, this means that the method is valid in the 
regime p < /i*, with /i* = 30/fc, 9/k 2 , 1.8/k 3 , 0.2/k 4 , 0.02/k 5 , 0.002/A; 6 for d = 4,5,6,7,8,9. 
Therefore, for k = 2 and d = 5 we get that our perturbative solutions describing two black 
holes on the cylinder are valid for p <C 2.2, in accordance with Fig. [TJ The values p* for k = 2 
black hole copies in 4 < d < 9 will be given in Table Q] in Section 16.11 

6 Further analysis of specific solutions 

In this section we analyze in more detail the two simplest multi-black hole configurations, 
namely two- and three-black hole solutions. This serves as an illustration of the general 
solution and its physical properties, but will also provide us with further insights into the 
structure of the phase diagram discussed in the previous section, including the possibility of 
existence of new lumpy black holes in Kaluza-Klein spaces. 

6.1 Two-black holes on the cylinder 

We start by examining the case of the two-black hole solution, i.e. we take a configuration of 
two black holes with mass fractions v\ = ^ + k and V2 = \ — k, where 0</t<l/2so that by 
convention M\ > Mi- Hence, k = corresponds to a configuration with two black holes of 
equal mass, while the limiting case k = 1/2 is the single black hole solution. The locations of 
the black holes are chosen as z\ = and the location of the second black hole is denoted as 
z\. For the equilibrium configuration we clearly have z\ = tt so that the two black holes are 
on opposite points on the circle. 
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We first focus on the equilibrium configuration. To compute the various thermodynamic 
quantities we need A^ 1 ' 2 ) defined in (|2.20p . which are given explicitly for the two-black hole 
case in Eq. (|A.9p . Furthermore, the expression for the sum Ya=i v i&^ is given in (jA.lOp . 



The curve fj4.9f) in the phase diagram is thus given by 

(d-2)Q(d-2 



n(ii;K) 



id-4 



(d-i)n d _i 



1 



An 2 1 1 



(6.1) 



Since the constant of proportionality is a monotonically increasing function of k one sees here 
explicitly that the inequality (|5.ip at k = 2 is obeyed, so that the slope in the n) phase 
diagram is bounded by that of a single black hole and two equal mass black holes. 
From (|4.1ip and (|A.9P we find the total entropy is given by 
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(6.3) 



where we used that ^(/u;^) = Si(fi]— k). In particular, we find from this the corrected 
entropy of one black hole on a circle <Sibh(m) = S(fi; 1/2) and that of two equal mass black 
holes S2eBu(^) = 0). We can now consider S(fi; k) for fixed total (rescaled) mass \i as n 
ranges from to 1/2. Physically, we expect that this is a monotonically increasing function of 
K since it should be entropically favored to have all the mass concentrated in one black hole, 
and the solution with two black holes is in an unstable equilibrium. As shown in Fig. [21 this 
is indeed the behavior we find when the mass of the system is not too large. 




Figure 2: Plot of the total entropy S of an equilibrium two-black hole configuration as a 
function of its mass distribution k, for a fixed total mass [i. This is a schematic plot for 
fj, < ij, c . 

We can in fact use the physical criterion that S(/j,; k) be a monotonically increasing function 
of k to get an upper bound fi c on the mass, in order for our perturbative approach to be valid. 
By examining the function (|6.2p in detail, we find that a condition that can be used to 
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determine the critical mass is 



d 2 S(jjL = fi c ;K) 



8k 2 



0. 



(6.4) 



K=0 



where above we also used that (8S(fJL, k) / dn)\ K= o = for all fi. Using the explicit expression 
()6.2p we have analyzed this equation for 4 < d < 9 and the results for \x c are listed in Table [TJ 
A necessary condition for our method to be valid is thus /i < fi c . We expect that for /i fi c 
our perturbative solution for the two-black hole configuration is valid. As illustrated in Table 
[TJ this is a less restrictive bound than the one found in the end of Section [5] based on a less 
precise consideration. 



d 


4 


5 


6 


7 


8 


9 


Pc 


14.4 


7.1 


3.9 


2.0 


0.97 


0.44 


p* 


15 


2.2 


0.2 


1 x 10" 2 


8 x 10" 4 


3 x 10~ 5 


fJ-GL 


3.52 


2.31 


1.74 


1.19 


0.79 


0.55 



Table 1: The upper bound /i c , imposed by entropy arguments, on the mass for the validity of 
the perturbative two-black hole results. For comparison the bound //* (see end of Section [5]) 
is shown along with the Gregory-Laflamme masses /iQL (see e.g. [12J). 



It is also useful to examine the temperatures of each of the black holes as we increase the 
mass. Clearly, for two black holes of unequal mass the zeroth order temperatures are different, 
and the system is not in thermal equilibrium. However, we can calculate the effect of the 
redshift on the ratio of temperatures, and examine whether this effect tends to equilibrate 
the black holes as we increase the total mass of the system. Using the first-order corrected 
temperatures (|4.2p and the expressions (|A.9P for A^ 1,2 ** one finds 

T 2 /1 + 2k\^\ AKC(d-2) , j * N 9 J , x 

Tr\T^) l 1 -t4=wr p ^" 1) " +0( ^J ■ (6 ' 5) 

where we eliminated po in favor of [i using (|4.7p . (|4.9p . 

For two unequal mass black holes (with M\ > M2) we have 0<k<1/2so the pre-factor 
in (|6.5p is greater than one. We now observe that as one increases the total mass \x the linear 
factor in fj, will be smaller than one, causing the ratio T2/T1 to decrease towards one. We thus 
conclude that the first order redshifts combine in such a way that increasing the total mass 
of two unequal mass black holes causes the temperatures of the two black holes to approach 
each other. This indicates that it may be possible that in the full non-perturbative regime 
the temperatures converge to a common value at the merger. 

Finally, we study the entropy of the more general configuration of two black holes without 
requiring the system to be in mechanical equilibrium. The total entropy S(/jl; k, z|) is obtained 
by using again (|4.1ip to compute the individual entropies Si^, but now substituting the z\- 
dependent functions A^ 1 ' 2 )^) given in Eq. (|A.8p . We consider then a fixed total mass [i 
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and mass distribution k, and vary the location z\ of the second black hole where < z\ < 
7r. Physically we expect that S([j,;k, z%) is a monotonically decreasing function of zj, with 
minimal entropy when the black holes are farthest apart and maximal entropy when they 
have merged into a single black hole. This is indeed the shown in Fig. [3j 

St 




Figure 3: Plot of the total entropy S of a two black hole configuration with fixed total mass 
\x and fixed mass distribution (here k = 0) as a function of the relative distance z% between 
the two black holes. We use a values of \i that lies below the critical mass \i c listed in Table 

m 

We can view the decrease of time evolution process in which two black holes 

initially separated by a distance tt on the circle are perturbed and then collapse into a single 
black hole. As seen in Fig. [3] the total entropy increases during this process, but the entropy 
diverges as the distance between the black holes goes to zero. This is expected since fields 
diverge when we let the distance between sources go to zero in the point-particle limit, and 
indeed A^' 2 \^) in (|A.8|) diverge as z\ — > 0. However, for physical sources, the minimum 
distance of approach between the sources is given by their size. In our case, a good estimate 
for this critical distance is given by the horizon size of a (d + l)-dimensional Schwarzschild 
(spherical) black hole with total mass //, given by 

On the other hand, we can compute the distance z c at which the entropy curve S(fi; k, z%) 
crosses the entropy Sibh(aO °f a single black hole configuration (see Fig. [3]), i.e. 

S(ji] k, Zc) = 5ibh(m) • ( 6 - 7 ) 

Comparison of the two critical distances p s and z c now provides an important check on the 
validity of our perturbative method, since we expect these two numbers to be of the same 
order. As illustrated in Table [2] this match indeed occurs with p s > z c , where for definiteness 
we have chosen k = 0. 

We thus conclude that also for non-equilibrium configurations the corrected thermody- 
namics leads to physically sensible results. 
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0.01 
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0.17 
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0.53 


1.13 


1.97 


0.54 


0.95 


1.68 


2.25 


Ps 


0.081 


0.26 


0.82 


2.58 


0.29 


0.62 


1.33 


2.54 


0.56 


1.00 


1.78 


2.34 



Table 2: Comparison of the two critical distances z c and p s in the case k = for some repre- 
sentative values of d and p (taken below the mass p c for which we can trust the perturbative 
results, see Table [1]). z c is the minimum distance imposed by entropic considerations, as il- 
lustrated in Fig. [3l and p s is the size of a (d + l)-dimensional Schwarzschild black hole with 
mass \i. 

6.2 Three-black holes on the cylinder 

In this subsection we discuss some features that can be addressed when we have three (or 
more) black holes, and we skip properties that are already present in the two-black hole 
configuration. In particular, by studying merges of two black holes we find evidence for new 
"lumpy" black hole configurations. 

For definiteness, take a symmetric three-black hole configuration in equilibrium, located 
at the points z\ = 0, z% = it — y, and z% = it + y. We also adjust the masses Mi = V{M such 
that i>2 = V3 = |(1 — i.e. black hole 2 and 3 have equal mass. We now want to increase 
the total mass of the system while maintaining equilibrium. The black holes will thus increase 
in size and fill more and more of the free space in between them. The question we want to 
address is whether the two black holes 2 and 3 with the same mass will merge first, before 
merging with black hole 1, or whether black hole 1 will merge with the other two before 2 
and 3 can merge. 

As in the previous subsection, our answer to this question is limited by the fact that our 
formulae are strictly valid only for small black holes interacting via Newtonian gravity, while 
the black hole merging process we wish to consider is certainly one where the full nonlinearities 
of Einstein's equations are important. However, we expect that with the available construction 
we can gain useful insights into the behavior of the system, so we proceed to examine this 
situation keeping in mind potential caveats. 

The question above can be addressed by analyzing the ratio 

X= Ps(1)+Ps(2) Z23 , (6.8) 

Z 12 Ps(2) + Ps(3) 

where p s ^ is the Schwarzschild radius of the i th black hole (defined as in (|6.6f) ) and zij is 
the distance between the i th and j th black holes. It is not difficult to see that this ratio is 
appropriate. Indeed, if black hole 1 joins 2 (and 3, by symmetry) first then at the point they 

merge one has = 1 and ^ > 1, so that X > 1. On the other hand, if 2 and 

3 merge first then at one has X < 1 at the merging point. 

We can express the ratio X defined in (|6.8p as a function of the distance y between black 
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hole 1 and 2 (and 3) as follows. First one uses the relation p B u\ oc {viM)~^ between the 
Schwarzschild radius and the black hole mass in d + 1 dimensions along with the fact that 
v 2 = ^3> so that 



Ps(l) + Ps(2) 
Ps(2) + Ps(3) 



1 + 



l 

d-2 



1 + 



^3 
^12 



l 

d-2 



(6.9) 



where we used the equilibrium conditions (|3.T[) in the last step. Finally, we substitute the 



explicit expressions (|3.4p for Vij where 212 



7r — y in V12 and 223 



2y 



in V23- Note that equilibrium requires Z23 < vr, so we only consider < y < tt/2. Equilibrium 
also demands that v-i = ^3 < ui/2. 

Collecting results, we use (|6.9|) to write the ratio in (|6.8p as 
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(27T) 



d-l 
d-2 



d-l 



(6.10) 



(2y)"- 1 [C (d-l,l-jf)-C(rf-l,l + if)] 



(2tt) 



d-l 



(tt-^-i [C (d- 1,1 -^)- C (d- 1,1 + ^)] 



We can understand (|6.10p as follows. A given value of y fixes the location and mass distribution 

of the system in equilibrium. Now let the total mass of the system increase. There is a critical 

value, call it M23, above which 2 and 3 are merged. Similarly above a critical value M12, 1 is 

1 

merged with 2 (and 3). X can then be expressed as the ratio (M23 / M12) d ~ 2 ■ Thus, if X < 1, 
as we increase the total mass of the distribution, black hole 2 and 3 will merge before 1 joins 
them, and vice-versa. 

A numerical analysis of (|6.1Up shows the following features, see Fig. [H For ^ < y < ^, 
one has X > 1; for < y < 5, X < 1; and for < y < y*, X > 1 again. Here, y* is a 
critical value that depends on the dimension of the spacetime. For example, for d = 4 one 
has y* ~ , while for d = 9 one has y* ~ . More generally, as d grows y* decreases and 
the interval where X < 1 grows. 




Figure 4: A typical plot of X — 1 versus the distance y ranging from to tt/2. 

To interpret these results first note that the case y = tt/2 describes a two-black hole 
configuration that is the limiting case of the three-black hole configuration where M\ — * 0, 
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and we have two equal black holes each with mass M/2 located at z = tt/2 and z = 3tt/2. The 
case y = it/3 corresponds to a symmetric configuration with three equal black holes equally 
spaced along the circle. The case y = yields essentially the single-black hole limit of the 
three-black hole configuration where M\ = M is centered at z = while M2 = M3 — > 0. 
Therefore, as y goes from to tt/2, the masses M2 = M3 increase from up to M/2, while 
M\ decreases from M to 0. 

Keeping these features in mind, and that § < 2/ < § implies Z23 > z\2, it follows that 
as the total mass increases black hole 2 and 3 will merge with black hole 1 before they meet 
each other. We thus expect, as observed above, that X > 1 for these values of y and X = 1 
at the boundaries of the range. To understand the behavior of X for < y < ^ it does not 
suffice to use purely geometrical arguments. Indeed, since Mi > M^ in this branch, black 
hole 1 seems to approach 2 (and 3) faster than 2 and 3 approach each other, but since at the 
same time Z23 < Z12, we should use the numerical analysis of X described above to determine 
what happens. This tells us that black hole 2 and 3 merge first, at least for y* < y < ^ where 
we found X < 1. However, for y values smaller than y* the numerical results for X are not 
reliable anymore, since in particular as y — > we see that X — > 00 which is due to the fact 
that our formulae are strictly valid in the point-particle limit where fields diverge when the 
distance between sources vanishes. 

The results above suggest that it could be possible that after the merging of the two black 
holes (2 and 3) we end up with a "lumpy" black hole (i.e. a 'peanut-like' shaped black object) 
together with an ellipsoidal black hole (1). It is conceivable that such a configuration would 
be a new static black hole solution in asymptotically A4 d x S 1 spacetimes. Generally if two 
black holes were to merge in this way, we expect that the resulting configuration would be 
singular. The singularity would arise if the surface gravities or temperatures of the two black 
holes differed, following standard results of \36\ I37|. In the above construction, however, we 
chose M2 = M3 to make the surface gravities identical in the merger. 

To discuss this further, note first of all that it is still true that the area of one spherical 

black hole of given mass is bigger than the sum of the areas of two isolated black holes with 

the same total mass. Nevertheless, the following argument suggests the possibility of lumpy 

objects for d > 4. In general dimension d+ 1, the horizon radius of a Schwarzschild black hole 

1 1 
scales as p s ~ M d ~ 2 , so starting from two black holes with /0 s (i,2) K M^ 2 we have at the 

1 1 

merging point a total radius p s ru + P s (2) x M^~ 2 + M 2 d ~ 2 . On the other hand a single black 

1 

hole with mass M\ + M2 has a radius p B n2) oc (Mi + M2) d ~ 2 ■ In four dimensions (d = 3) 
this scales the same way as the total radius of the merged object so we expect the formation 
of a spherical black hole [38J. However, for d > 4 the power in the exponent is less than one 
so that p s n) + p s (2) > Ps(i2) • Hence we should expect that the resulting merger configuration 
will not be surrounded by a spherical horizon, as would occur for d = 3. 

As a consequence, it seems that for d > 4 this analysis does not rule out the possibility 
of having a configuration of a lumpy black object with "centers" at z = it ± y kept in an 
unstable equilibrium by a black hole at z = (and the respective copies). Note also that 
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the argument above suggests that the higher the dimension, the more likely it is that lumpy 
black holes will exist. Finally, we emphasize that the above analysis should be read within 
the earlier-mentioned limitations of our construction. 

Note that in asymptotically flat space new stationary black holes with similar 'rippled' 
horizons of spherical topology have been argued to exist in Ref. [39] by considering ultraspin- 
ning Myers-Perry black holes in dimensions greater than six. While in that case the ripples 
are supported by the angular momentum J in our case they are supported by the external 
stress of the other (ellipsoidal) black hole. It would be interesting to generalize the analysis 
above to configurations with more bumps, for example taking a symmetric four-black hole 
configuration with M2 = M3 = M4. 

7 Conclusions and outlook 
7.1 Summary 

In this paper we constructed solutions of the vacuum Einstein equations describing multi- 
black hole configurations on the cylinder x S 1 with d > 4, in the limit of small total 
mass, or, equivalently, in the limit of a large cylinder. These solutions generalize the analytic 
solutions found for the single black hole on the cylinder \VI \ 118 1 H9 l l20 | l21j. Furthermore, they 
generalize the so-called copies of the single black hole solutions corresponding to having equal 
mass black holes distributed equidistantly around the cylinder [281 E2] • The new solutions are 
valid to first order in the total mass, and are constructed using the technique of [T7] based on 
an ansatz for the metric found in |16j . 

Using the first-order corrected metrics for the multi-black hole configurations we have 
studied their thermodynamics. Included in this is one of the central results of this paper: The 
relative tension (binding energy) n as a function of the total (rescaled) mass fj,, as given by 
Eq. (|4,9p . Using this, we have shown how the solutions appear in the (/i, n) phase diagram 
[12[ I29j. together with the other known solutions that asymptote to A4 d x S 1 . We observed 
that a multi-black hole configuration with k black holes has k independent parameters. This 
implies a continuous non-uniqueness in the (jJ>,n) phase diagram (or for a given mass), much 
like the one observed for bubble-black hole sequences [30] . 

The multi-black hole configurations have to be in mechanical equilibrium in order to have 
a static solution. We have identified where this requirement appears in the construction of the 
solution, and we have furthermore examined how to build such equilibrium configurations. 
Moreover, we have described a general copying mechanism that enables us to build new 
equilibrium configurations by copying any given equilibrium configuration a number of times 
around the cylinder. 

Finally, we examined in detail configurations with two and three black holes. For two 
black holes we verified the expectation that one maximizes the entropy by transferring all the 
mass to one of the black holes, and also that if the two black holes are not in mechanical 
equilibrium then the entropy is increasing as the black holes become closer to each other. 
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These two facts are both in accordance with the general argument that the multi-black hole 
configurations are in an unstable equilibrium and generic perturbations of one of the positions 
will result in that all the black holes merge together in a single black hole on the cylinder. For 
the three black hole solution we examined and found preliminary evidence for the hypothesis 
that for certain three-black hole configurations two of the black holes can merge into a lumpy 
black hole, where the non-uniformities are supported by the gravitational stresses imposed by 
an external field. 

From the first-order corrected temperatures one can show that the multi-black hole config- 
uration are in general not in thermal equilibrium. The only configurations that are in thermal 
equilibrium to this order are the copies of the single-black hole solution studied previously 
[28} [29l [T7] . As a further comment we note that Hawking radiation will seed the mechanical 
instabilities of the multi-black hole configurations. The reason for this is that in a generic con- 
figuration the black holes have different rates of energy loss and hence the mass ratios required 
for mechanical equilibrium are not maintained. This happens even in special configurations, 
e.g. when the temperatures are equal, because the thermal radiation is only statistically uni- 
form. Hence asymmetries in the real time emission process will introduce disturbances driving 
these special configurations away from their equilibrium positions. 

7.2 Discussion of the phase structure 

We now examine the appearance of our new multi-black hole phases in connection to the known 
phases of black holes and black strings on the cylinder (See [HEH3] f° r reviews). In particular, 
as mentioned in the introduction there is the well-known phase of the uniform black string 
(UBS) as well as the non-uniform black string (NUBS), emanating from the uniform phase at 
the Gregory-Laflamme point fJ-GL- Recently, numerical investigations [TT | \13 \ \T Q \TE[ [22l [23j 123] 
confirmed the prediction [25] that the non-uniform phase connects via a horizon topology 
changing phase transition [251 E3 US] to the phase of a single localized black hole (LBH) 
(see Fig. [TJ. This point is generally referred to as the merger point. 

Moreover, as reviewed in Section[3l for any solution that falls into the SO(d— l)-symmetric 
ansatz (|2.ip of Ref. [16], one can obtain a copied solution [281 ESI Q2] by changing the peri- 
odicity of the circle from L to kL with k an integer. As mentioned above, this includes the 
localized black hole phase, from which one generates in this way the multi-black hole solutions 
with k equal mass black holes, which we denote by LBH^. It also includes the non-uniform 
black string phase, from which we generate copies which we denote by NUBSfc, emerging from 
the uniform phase at critical mass /xgl/^ -3 - This thus means that the LBH/% phase will 
connect to the NUBS^ phase via a horizon topology changing phase transition at the fc-copied 
merger point (see Fig. [I] for k = 2). 

We now turn to the question posed in Section [5} Where do all the new multi-black hole 
phases end in the phase diagram? For definiteness, let us consider again configurations with 
two black holes. The LBH and NUBS phases are connected via the topology-changing merger 
point, and likewise the LBH2 and NUBS2 phases are connected via the 2-copied merger point. 
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As explained in Section [5] all two-black hole configurations with unequal mass lie (at least for 
small masses) in between these two limiting phases and it is not clear where these phases will 
end up in the phase diagram. Two scenarios where given in Section [5j and we now examine 
in more detail the possibility of the second scenario, namely that the black holes merge into 
a new non-uniform string. Recall that this would require the temperatures of the black holes 
to approach each other at the merger point. 

First of all, we have seen in Section 16.11 that our first order result for the temperatures 
shows that the temperatures of the two black holes are redshifted in such a way that they 
tend to approach each other. This lends credibility to the possibility that indeed in the full 
non-perturbative regime the temperatures may converge to a common value at the merger. 
If this is the case, it seems to suggest that there would exist new non-uniform black strings 
beyond the NUBSfc phases, to which the unequal mass black hole configurations could connect 
via new merger points 

If smooth mergers do occur for different size black hole configurations, an important 
question to consider is whether this a generic feature, or if it only happens for particular 
configurations. Consider for the example the case of two black holes, for which we have two 
free parameters namely the total mass and the ratio of the individual masses. Equating their 
temperatures fixes the mass ratio as a function of the total mass (see Eq. (|6.5|) ) . Similarly 
achieving a merger of the black holes also fixes the ratio with another function of the total 
mass. If we imagine these two functions of the mass to be independent, it follows that we only 
expect these two functions to intersect at discrete points in the space of parameters defining 
the configuration. On the other hand, if these two functions are not independent, due to the 
interrelation between geometry and energy in General Relativity, one can instead imagine 
that the two functions always intersect, so that the smooth mergers are a generic feature. 

As discussed above, smooth mergers for different size black hole configurations suggest that 
new non-uniform black string phases exist. If this is the case, there are certain constraints 
on such new phases from general arguments. Firstly, it is clearly not possible that there are 
non-uniform black strings emerging from the uniform black string in the range /iGi,/2 d ~ 3 < 
A* < ^GL [S]- Also, it does not seem possible that one can have other branches than the known 
ones coming out of the Gregory-Laflamme point (or its /c-copies) of the uniform black string 
given the higher-order perturbative analysis of Ref. [TO]. Secondly, it is impossible to (locally) 
have a continuum of non-uniform black string solutions in the phase diagram. To prove this 
assertion imagine that there is a two-dimensional continuous parameter space of solutions and 
consider two points, say A and B, in this continuum. It follows from the continuity that one 
can always connect these two points by two different paths of solutions. Imagine now that 
the two-dimensional space of solutions projects into a two-dimensional region in the (//, n) 
phase diagram. If we then furthermore take the paths so that n in path 1 is greater than n 
in path 2, then we get an contradiction when using the Intersection Rule of Ref. |12j . This is 

10 Note that the original argument by Kol [251 [T] for the merger transition of the LBH and NUBS phases was 
based on Morse theory, which loosely speaking implies that the LBH phase cannot end in "nothing" . 
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because c^Si/S^) = — n^MbMj ((d — l)TiT2Sf) where the indices on the quantities refer 
to the paths. Since n\ > 712, the right hand side is strictly positive. Thus, the ratio S1/S2 in 
point B should be greater than 1, but that is not possible since the two paths should go to the 
same solution. We thus conclude that a locally continuous space of solutions is impossiblj^l. 
except in the very special case where the continuous space of solutions projects onto a one- 
dimensional subspace in the (fi, n) phase diagram. This provides a further argument that 
smooth mergers would only occur at discrete points, because there could only be a discrete 
set of non-uniform string solutions to which the the merging black holes could connect. 

Given these two constraints, there is still the possibility that new non- uniform black strings 
may exist. Namely, it is conceivable that the NUBS^ phases (k > 1) develop their own zero 
modes as one moves away a finite distance away from the GL point (or its A;-copies). This 
is a non-perturbative effect that would not show up in the perturbative analysis of Gubser. 
These zero-modes on the non-uniform black string would in fact imply that they have some 
region in which they are respectively classically unstable or stable, just as for the uniform 
black string. Such a bifurcation of new non-uniform strings from the presently known ones 
would also be discrete and thereby evade the second restriction presented above. Furthermore, 
in this scenario one could imagine a fractal structure of further bifurcations into new non- 
uniform strings, all of which eventually end up in a particular multi-black hole configuration. 
If true, this would fit well with the smooth mergers of different size black hole configurations 
occurring only at a discrete points in the space of configurations. It would be very interesting 
to explore this possibility further. 

Another point that we already alluded to in Section [5] is that we do not expect the phases 
of two black hole configurations to stay in between the LBH and LBH2 curves in the (/u, n) 
phase diagram (and similarly for multi-black hole configurations with more than two black 
holes). To see this consider the LBH curve in Fig. [TJ This curve has a point at which /i is 
maximal, occurring well before the merger point. Beginning with this maximal mass single- 
black hole configuration we can add a tiny black hole on the opposite side of the circle and 
reach a two-black hole configuration with greater mass than the original configuration. This 
clearly implies that the two-black hole configurations can extend outside the wedge bounded 
by the LBH and LBH2 curves. In fact, one can similarly argue by starting from the extremal 
point on the LBH2 curve that by removing a tiny mass from one of the two black holes, 
one can reach a two-black hole configuration to the left of this curve. Another interesting 
example comes from adding a pair of tiny black holes to any LBH2 configuration to produce 
a four-black hole configuration in its neighborhood in the phase diagram, very far away from 
the wedge enclosed by the LBH3 and LBH4 curves. The above reasoning can be extended 
by imagining further additions of tiny masses, in more complicated starting configurations, 
leading to a intricate pattern of crossings of lines in the (fi, n) phase diagram. 

11 Note that implicit in the above argument is the assumption that there is only one connected horizon with 
a given temperature. Thus, the fact that multi-black hole configurations cover a continuous region in the phase 
diagram is not a contradiction because they contain disconnected horizons typically at different temperatures. 
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We have also presented evidence in this paper for the possibility of a new class of static 
lumpy black holes in Kaluza-Klein space. Again, it would be interesting to study this further, 
and examine how these in turn might connect to new non-uniform phases. 

7.3 A fluid analogy 

It is also interesting, though more speculative, to consider the appearance of the multi-black 
hole configurations in relation to an analogue model for the Gregory-Laflamme (GL) instabil- 
ity, recently proposed in Ref. [30]. There it was pointed out that the GL instability of a black 
string has a natural analogue description in terms of the Rayleigh-Plateau (RP) instability of 
a fluid cylinder. It turns out that many known properties of the gravitational instability have 
an analogous manifestation in the fluid model. These include the behavior of threshold mode 
with d, dispersion relations, the existence of critical dimensions and the initial stages of the 
time evolution^ ( see Refs. gDJ B21 g3] for details). 

Since our reasoning below relies on the time evolution of the system and its endpoint, it 
is worth mentioning that the full time evolution of the RP instability is well known (both 
numerically and experimentally, see Refs. |44[ 1451 1^61 1^1 f° r details). On the gravity side only 
the initial stages of the GL instability has been numerically studied so far [48J. Comparing 
with the fluid system there is an interesting match between the initial stage of the evolution 
in the two systems. Starting from a single sinusoidal perturbation both develop an almost 
cylindrical thread or neck in between the two half rounded boundary regions. This can be 
confirmed by comparing Fig. 1 of [16] (which describes the full RP evolution) and Fig. 2 of 
[4"5] (that describes the initial stage of the GL evolution). 

One should be cautious when applying the analogue model, especially in what concerns 
the evolution of the systems. The reason is that the analogy is partly based on the similarity 
between the first law of black hole thermodynamics and the fluid relation dE = TdA where 
E is the potential energy associated with surface tension (free energy), T the effective surface 
tension and A the surface area of the fluid. This means that both systems tend to extremize 
the area. However, on the gravity side we know that a black object evolves such that its 
horizon area never decreases, whereas a fluid evolves toward a configuration with smaller 
area, since this decreases its potential energy. Despite these reversed dynamical features, it 
is worthwhile to notice that just like a multi-black hole system will maximize its entropy by 
merging into one single black hole containing all the mass, so will an array of fluid droplets 
merge into a single drop in order to minimize its surface area at fixed volume. 

Having alerted the reader to these caveats, we proceed with the analogy in hand, consid- 
ering the time evolution of the fluid in further detail. A representative study of particular 
interest for our purposes was carried out in [16] . The main conclusion is that if we start with 

12 Recently, another feature of these instabilities has been matched. If rotation is added to the fluid the 
strength of the fluid instability increases because the centrifugal force is bigger in a crest than in a trough 
of the configuration. On the gravity side it was found in Ref. [4T] that rotating black strings, even for large 
rotation, are indeed also unstable to the GL instability. 
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a single sinusoidal perturbation in a cylindrical liquid bridge, the higher harmonics generated 
by non-linear effects are responsible for the development of a long neck that break^f] in a 
self-similar process [46 1 |4"4"1 145] . We end up with an array of satellite drops with different 
sizes. Hence, if the correspondence indeed extends to the full evolution, the multi-black holes 
would be the natural gravity analogues of the main drop and satellite droplets array observed 
in the fluid analysis. 

Furthermore the analogue model would thus argue in favor of the scenario in which the 
neutral black string will pinch off. Moreover, the multi-black hole configurations constructed 
in this paper would play an important role in the intermediate stages of the GL instabil- 
ity. It would be interesting to examine this application of the analogue fluid model and its 
consequences more closely. 

7.4 Outlook 

The study of Kaluza-Klein black holes and their high degree of non-uniqueness can be viewed 
in the broader context of studying black objects in higher dimensional gravity. Here, research 
in the last years has revealed that also in asymptotically flat space a very rich phase structure 
of stationary black objects is expected. In particular, in five-dimensional Einstein gravity there 
exists, beyond the rotating Myers-Perry black hole, a black ring solution (see |5L)j for a 
review). Recently further new stationary solutions, called 'black Saturns' [32l [33] , have been 
constructed explicitly in five-dimensional gravity. These solutions, consisting of a spherical 
black hole with black rings around it, are similar to the multi-black hole configurations, in 
that the generic solution is not in thermal equilibrium, with different temperatures for each 
connected component of the event horizon. Furthermore, one may compare the configurations 
with highest entropy in the two systems. It was shown in Ref. [M] that the maximal entropy 
configuration for fixed mass and angular momentum consists of a central, close to static, black 
hole and a very thin black ring around it. For any value of the angular momentum, the upper 
bound on the entropy is then equal to the entropy of a static black hole of the same total 
mass. These maximal entropy black Saturns are not in thermal equilibrium. In some sense 
the same features are observed for multi-black hole configurations. If we restrict to the case of 
two black holes, the highest entropy configuration (see Section |6.1|) is that of an infinitesimally 
small black hole together with a large black hole, i.e. far away from thermal equilibrium. The 
entropy of that configuration is bounded from above by that of a single black hole of the same 
mass. 

It is also worth emphasizing that the solution technique employed in this paper can be 
applied to other black hole systems where one lacks the symmetries or other insights to 

13 Note that at the pinch-off there is another similarity that characterizes both instabilities. On the gravity 
side, one would need to use quantum gravity when the pinch-off region reaches the Planck scale and General 
Relativity is no longer valid. Likewise, close to breakup of the fluid, when the radius of the liquid bridge is of 
molecular size, the (continuum) hydrodynamic theory is no longer a good approximation and simulations of 
the molecular dynamics are required. 
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construct exact solutions. The general idea is to identify a suitable perturbation parameter 
of the putative solution, and follow similar steps as outlined in Section [2 

Another open direction to pursue is to apply numerical techniques to extend the construc- 
tion of multi-black hole configurations into the non-perturbative regime, as was successfully 
done for a single black hole on a cylinder in five and six dimensions \22\ [23j 124j . Such an 
analysis could confirm whether indeed there are multi-black hole configurations for which the 
temperatures converge when approaching the merger points as one increases the mass, as was 
discussed in Section 16.11 Furthermore, it is possible that in this way one could confirm the 
existence of the lumpy black holes conjectured in Section HOI where we recall that these are 
most likely for higher dimensions. 

A further, but technically complicated, direction to pursue is to extend the solutions of this 
paper to the next, i.e. second, order. For the case of a single black hole in five dimensions, the 
second order correction to the metric and thermodynamics have been studied in [20]. More 
generally, the second order correction to the thermodynamics was obtained in Ref. [21] for all 
d using an effective field theory formalism in which the structure of the black hole is encoded 
in the coefficients of operators in an effective worldline Lagrangian. It would be interesting to 
obtain the second-order corrected metric and thermodynamics for the multi-black hole case 
considered in this paper. 

There are also potential applications related to string theory and gauge theory. It is known 
that the phases of Kaluza-Klein black holes are related via a boost /U-duality map [51] (see 
also [52, 53j) to phases of non- and near-extremal branes on a transverse circle, appearing 
as solutions in type II string theory or M-theory. Via the gauge/gravity correspondence 
[54\ [55] this has implications for the phase structure of the dual non-gravitational theories at 
finite temperature. For instance, it is possible to obtain in this way non-trivial predictions 
|53i I5TI I5S] about the strong coupling dynamics of supersymmetric Yang-Mills theories on 
compact spaces and of the thermal behavior of little string theory. 

As an important example, Ref. [53] considered finite temperature two-dimensional super- 
symmetric Yang-Mills on a spatial circle, which by the boost/U-duality map is related to 
the phase structure of Kaluza-Klein black holes in ten dimensions. The corresponding phase 
structure that is present at strong coupling in the two-dimensional Yang-Mills theory on the 
torus Sp x S 1 was then qualitatively matched to the phase structure in the weakly coupled 
gauge theory. In particular, it was found in [53] that the eigenvalue distribution of the spatial 
Wilson loop distinguishes between the three different phases seen at strong coupling: The 
uniform phase corresponds to a uniform eigenvalue distribution, the non- uniform phase cor- 
responds to a non-uniform eigenvalue distribution and the localized phase maps to a gapped 
eigenvalue distribution. It would be interesting to see if there are also multiply gapped eigen- 
value distributions (see e.g. Ref. [51]), corresponding to the localized phase of multi-black 
holes found in this paper. While those would probably be unstable as mentioned above, they 
may still appear as unstable saddle points. 

Finally, we remark on an open direction that is related to microscopic calculations of the 
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entropy of black holes. In Ref. [58] (see [59] for a short summary) the boost/U-duality map of 
|51j was extended to the case of branes with more than one charge. One of the results is that 
by starting with neutral Kaluza-Klein black holes in five dimensions one can generate five- 
dimensional three-charge black holes on a circle, obtained from corresponding three-charge 
brane configurations in type II/M- theory via compactification. In particular, when one applies 
this map to a single neutral localized black hole one obtains a three-charge black hole localized 
on the transverse circle. For this case, it was shown that in a partial extremal limit with two 
charges sent to infinity and one finite, the first correction to the finite entropy is in agreement 
with the microscopic entropvF^I By applying the map to the multi-black hole solutions of this 
paper one will generate three-charge multi-black holes on a circle. The results of Section H] 
can then be used to compute the first correction to the finite entropy of these three-charge 
multi-black hole configurations, and it would be interesting to then derive these expressions 
from a microscopic calculation as well. Furthermore, in Ref. [60J a simple microscopic model 
was proposed that reproduces most of the features of the phase diagram of three-charge black 
holes on a circle, including the new non-uniform phase. It would be interesting to see if this 
model can also account for the corresponding localized three-charge multi-black hole solutions. 
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A Data for two unequal mass black holes 



In this appendix we give some useful explicit expressions for the quantities that are involved 
in the construction of the two-black hole configuration, further discussed in Section 16.11 

The mass fractions of the two black holes are taken as v\ = | + k, = \ — k and the 
equilibrium configuration is chosen such that the first black hole at z* = and the second 
at z\ = 7r, i.e. at opposite points on the circle. The function (|2.7|) entering the Newtonian 
potential is then given by 

F(r, z) = 2 d - 3 F(2r, 2z) + 2 d ~ 2 K T{2r, 2z) , (A.l) 

where we have defined 

Hr,z)= " O"' Hr,z)= — O"- ( A " 2 ) 

m =_oo [r 2 + {z- 2vrm) 2 ] 2 m =-oo [r 2 + (z - 2vrm) 2 ] 2 

The function J-(r, z) is identical to the one entering the Newtonian potential for the single 
black hole case, and details can be found in Appendix B of [16]. Using Poisson resummation 
the large r expansions of the two functions in (|A.2p are obtained as 

T(r, z) = rf d 3 (l + 2 Y^ h{mr) cos(mz) J , (A. 3) 

m=l 



2fe x 



H r > z ) = ~d=z Yj h ^ m + V 2 ] 7- ) cos([m + 1/2]*) , (A.4) 



r 

m=0 

where and /i(x) are defined in (|2.9p . (|2.10p respectively. 

Note that for even d, the relevant Bessel function takes the form of a polynomial of finite 
degree. This allows to write explicit expressions for J-(r, z) and J-(r, z) (and similarly for 
v(r, z) in f)2 . 14|) ) . For example, for d = 4 one has 

. . 1 sinhr A , , 1 sinh(r/2) cos(2;/2) „ . , . 

F(r,z) = - : , T(r,z) = \L2 K -U- , ford = 4. A.5 

2r cosn r — cos z r cosn r — cos z 

For d = 6 one has 

1 sinhr + 1 sinh 2 (r/2) cos 2 (z/2) - cosh 2 (r/2) sin 2 (z/2) 

4r 3 cosh r — cos z 2r 2 (cosh r — cos z) 2 

p 1 sinh(r/2) cos(z/2) ^ 1 sinh 2 (r/4) cos 2 (z/4) - cosh 2 (r/4) sin 2 (z/4) 

2r 3 coshr — cosz 8r 2 [cosh(r/2) — cos(z/2)] 2 

1 cosh 2 (r/4) cos 2 (z/4) - sinh 2 (r/4) sin 2 (z/4) „ . A , 

+ 7^-^- f n — t\ ; — ; — TTo i for d = 6 . (A. 6 

8r 2 [cosh(r/2) + cos(^/2)] 2 ' v ' 

For the small p expansion, we first present the results for general location z\ of the second 
black hole, restricting to the equilibrium configuration = 7r at the end. In the region near 
the first or second black hole respectively, we have from (|2.19p the expansions 

F ^ z ) - y=r + A(1) > F ^ z ) - y^r + A(2) > (A. 7) 
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where Ad' 2 ) are computed from (|2.20p 
Ad) 



x) _ (l + K )2({d-2) I-« ( \-k 



d-2 



(27T) 

2) _ (i-«)2C(d-2) 



(4)^2 (27T) 



d-2 



(2JT) 



d-2 



+ 



(27T - 



*W-2 



+ 



i d-2 



(2vr 



2tt 



Cjd-2,2-^- J+Cld-2,^- 



2vr 



(Ai 



and we recall the definitions (|2.23j) . (|2.24p . In particular, for the equilibrium configuration 
z% = 7T these expressions reduce to 

A (D = [2 d-3 + 2k{1 _ 2 d- 3)] gCg_l2) > A (2) = [2 d-3 _ 2k{1 _ 2 d- 3)] 2C(d ~ 2) (A _ g) 

Note that for k = 1/2, the expression for Ad) reduces to the correct result for a single black 
hole. Finally, we record the sum 



(27T) 



d-2 



l + 2 5 - d 1 -2 d - 3 ) K 2 



(A.IO) 



which is used in the text to compute various thermodynamic quantities. 
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